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In turbulent flows the n'th order structure functions S„{R) 
scale like when R is in the "inertial range". Extended 
Self-Similarity refers to the substantial increase in the range 
of power law behaviour of S„ (R) when they are plotted as a 
function of 52 (R) or 5*3 (-R) . In this Letter we demonstrate this 
phenomenon analytically in the context of the "multiscaling" 
turbulent advection of a passive scalar. This model gives rise 
to a series of differential equations for the structure functions 
Sn{R) which can be solved and shown to exhibit extended self 
similarity. The phenomenon is understood by comparing the 
equations for S„{R) to those for S„{S2)- 

PACS numbers 47.27.Gs, 47.27.Jv, 05.40. +j 



The fundamental theory of turbulence concerns itself 
with the universality of the small scale structure of turbu- 
lent flows and their characterization by universal scaling 
laws ||]. Since the days of Kolmogorov ||^ one expects 
the structure functions to scale with R as long as 

R is in the inertial range 



for ri<^R-^L 



(1) 



Here 77 and L are the viscous and outer lenghts re- 
spectively. For a given turbulent field u(r, t) the 
n'th order structure function is defined as Sn{R) = 
{[u{y -f- R, i) — M(r, t)]") where the angular brackets (. . .) 
denote an average over r and t. Komogorov's 1941 pic- 
ture of turbulence (K41) predicted that the exponents 
of the structure functions of the velocity field satisfy 
Cn = n/Z. Further research along the same lines sug- 
gested that the structure functions of passive scalars den- 
sities also have the same exponents On the other 
hand experiments 1^,^ have indicated deviations from 
this scaling and even the possibility of multiscaling, i.e. 
a non-linear dependence of C„ on n. These findings raise 
fundamental questions as to our understanding of the 
universality of the small scale structure of turbulence. 

Unfortunately, the experimental measurement of the 
scaling exponents in turbulent flows is hampered by the 
fact that the range of scales that exhibit scaling be- 
haviour is rather small for Reynolds numbers Re acces- 
sible to date. Even though the viscous scale rj decreases 
like Re~^/^, the scaling behaviour (|l|) is expected to set 
in only at about IO77, and cease at about L/100. Thus 
even at Re ^ 10^ the inertial range is not broad enough 
to allow accurate measurements of Cn, particularly for 
larger values of n. 



A partial resolution of these difficulties has been re- 
cently offered by changing the way that the experi- 
mental data is presented. Instead of plotting log-log plots 
of Sn{R) vs. i?, it was argued that plotting log- log plots 
of Sn{R) vs. S2{R) or S3{R) reveals much longer scafing 
ranges that can be used to fit accurate exponents. This 
method was dubbed Extended Self Similarity. While the 
effectiveness of this method has been demonstrated 
the reason for its success remains unclear. Therefore its 
reliability is not guaranteed. The aim of this Letter is to 
address this method theoretically for Kraichan's model 
of a passive scalar advected by a very rapidly varying 
velocity field | ^. In this model multiscaling is known 
to occur [|l^, 

and we have sufficent analytic under- 
standing of the mechanism of (multi)scaling to assess why 
Extended Self Similarity could hold. One should stress 
that in the case of normal scaling Extended Self Simi- 
larity can be rationalized rather straightforwardly [Q. 
To our knowledge no analytic derivation of this effect is 
available in a case exhibiting multiscaling. 

In this model the equation for the passive scalar con- 
centration T(x) is given by 



dT 
'dt 



(2) 



where the velocity field u is taken as Gaussian random 
and with a (5-function correlation in time. The structure 
functions S2n{R) are know to obey differential equations 

R'-''-^R''-'HR)-^S2n{R) = J2niR) ■ (3) 

Here h{R) = H{R/C)^>' for R< C where H is a. dimen- 
tional constant and C is some characteristic outer scale 
of the driving velocity field. The Scaling exponent C,h is a 
chosen parameter in this model, taking any desired value 
in the interval (0,2). The function J2n{R) is known ex- 
actly for R in the inertial range 0, but for our purposes 
we need to know it also in the near viscous range. A form 
that is appropriate in both ranges has been suggested by 
Kraichnan [l0|] : 

J2„(i?) = 2nn^l^[S/'S2{R) ~ V'S^m . (4) 

This form has been recently supported by numerical sim- 
ulations ||l^ . In the inertial range the last two equations 
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have a multiscaling solution for S2n{R) with the expo- 
nents 



C2n = i V(3-C2)2 + 12nC2 + C2-3 



(5) 



Our first step in studying the issue of Extended Self- 
Similarity using Eqs. (|3|j^) is to non-dimensionalize them. 
Rescale lengths by = C{ii/ HY^^'^ and the scalar 
density hy Q = ej k. The dimensionless functions, 
D2niy) are: 



They satisfy the equations 

d'^D2n , 2 + CridD2n 2nD2 



y = R/rd 



dy 



D2y^'' 



<Py 



(6) 



(7) 



The success of extended self similarity for this model can 
be assessed from Fig.l which compares plots of D2n vs. 
y with plots of D2n vs. D2. The data is based on the nu- 
merical solution of Eq.(|^) using the IMSL/IDL's Adams- 
Gear ODE-solver. The efficacy of Extended Self Similar- 
ity speaks for itself. 

Next we exhibit the phenomenon analytically by find- 
ing the first order dissipative correction to the inertial 
range scaling of the structure functions. In the inertial 
range the exponent of D2(y) is ^2 = 2 — We can 

find the first correction to the power law, going into the 
dissipation range, by writing 1)2(2/) — 2y'»^/3C2 -I- e(y)- 
We then find the following ODE for e{y): 



d^y 



y 



Cr,de ^ 4(1 



■Cn) 



dy Sy^^i 







(8) 



There are two possibilities for the leading order solution 
of this equation. If all the terms are of the same order 
the solution is ei = — 2y^~^''''/3(l -I- ^2)- If the first two 
terms are larger we have £2 = B2 + ei, where now ei 
appears as the next order correction after a constant B2 ■ 
One can see that the dominant solution for small is ei , 
with a transition to £2 at = 1- Thus the near- inertial 
solution for D2 is 



D2iy) 



2y^- 



3C2 



1 



C2y 



C2-2 



D2{y) 



2y^ 

3C2 



I + C2 
, C2 



7C2 



for C) < 1 , 
for Cn> I , 



(9) 
(10) 



where C2 is a constant. Note that for the case (^h = 2/3 
the ODE for D2{y) can be solved exactly in terms of 
elementary operations and the results are consistent with 
this solution. 

Next we repeat this procedure for D2n{y)- We substi- 
tute D2n{y) = ^2n?/'^""[l + ^iv)] ^ Eq.(|). For the case 
C/i > 1 we find 



^ + -^(2(2 

dy2 y dy ^ 



C2) 



3nC2C2 



y 



, (11) 



with a leading order scaling solution ^(y) = i?2nJ/~^^- 
The constant is found from Eq.(|ll|) and the final result 
for C2 < 1 is 



D2n{y) = A2„y'^^" 



1 



3nC2 



(2C2n 



2C2)2/f= 



(12) 



At this point we can assess the efficacy of Extended Self 
Similarity for the case C/i > 1. If all functions D2n were 
an exact power law in argument D2 we would have 



D2n (y) 



or in the same order 



2{y^' + C2 



3C2 



C2„/C2 



(13) 



D2n (y) = A2ny'''" 1 



C2nC2 

C,2y^^ 



(14) 



The ratio ri between the coefficients of y '■^ in Eqs. (^1 



and (14) measures the effectiveness of this presentation. 
Using Eq. (|^) one finds 



3nC2 



"C2 



C2n(2C2« + 3-2C2) 2nC2-C2„ 



(15) 



If ri = 1 we would have perfect Extended Self Similarity; 
ri = is what we would get if we just plotted S2n as a 
function of r. In Fig. 2(a) we can see a graph of ri as 
a function of C2 for n = 2, rt = 4 and n = 6. We can 
see that r is always in the interval [1/2, 1]: from Eq. (|^) 
lini„^oo ''I = 5, and from Eqs. (|5|),(p^ lim^^^o?'! = 1- 
The convergence to 1/2 is rather slow and not uniform 
in (.2- 

The main point of the Extended Self Similarity is that 
52 (which is a function of R) is a "better" variable than 
the distance R itself when it comes to scaling behaviour. 
Thus, to complete our analysis we want to write the equa- 
tions for S2n as ode's with respect to argument S2 with 
the inertial range scaling plus a remainder term. Chang- 
ing variables and writing the S2n as a function of 6*2, 
we expect to find that the remainder term, evaluated 
in terms of the inertial range estimates, is smaller than 
when using r as a variable. To this aim we introduce a 
new set of function F2n{x) of argument x — D2{y)'. 

F2n{x)=D2n{y), X = D2{y) . (16) 
Then Eqs. (0) gives for F2n{x) the following ODE 



d^F2n 



[/(X)]2-C2 



dF2n nF2n{x) 



dx 



(4 - C2) df 
f{x) dx 



= 



(17) 



where f{x) is the inverse function of 5*2(2/): f{x) = y if 
X = 82(1/) ■ Now let us introduce a third set of dimen- 
sionless structure functions, G2niz) such that 



2 



(18) 



Note that the functions F2n and G2n are identical in the 
inertial range, since one is a function of 5*2 and the other 
is a function of its power law in the inertial range. Indeed, 
it will be seen below that they satisfy the same ODE in 
the inertial range with a power solution x''^"/''^. Using 
(0) one derives the following ODE for G2n{x): 



2n 



3 

C2X [ dx 
3nG2nix) 



dG2n nG2n{x) 



(C2 



9{x) 

dx dx'^ 



= (19) 



with the function g(x) satisfying g(x = D2{y)) = D2{y). 

Consider first the C2 < 1 case, and use the near-inertial 
approximation for 5*2(2/), Eq. (0)- Inverting this relation 
we find the relevant approximation for / in terms of x, 



fix) 



^ K2X 



I/C2 



2C2 

3(1 X 



(20) 



Using this and introducing constants 



■'2n — 

Cn„x^^"^''^ in the remainder term we find that the near- 



inertial version of Eq. (uTY) is 



d^F2n 

d^x 



3 

C2X 



dF2„ nF2n [x) 



dx 



(21) 



2C2C2n 

73 IC2« 

^2 



nC,2)x 



C2„/C2-3 



Doing the same to Eq. (n9|) we have 



D^G2 

d^x 



dG2n _ nG2n{x) 
dx X 



3 

C2X 

2C2G2 



(22) 



CI 



J (2 







Indeed, we see that the first lines of these equations , 
which are the inertial range contributions, are the same. 
Also the second lines are of the same order in y but the 
ratio r2 between the coefficients in these lines (reminder 
terms) is 



r2 



C2n - 11(2 
71(2 



(23) 



Since the Holder inequalities imply that C,2n < 1^(2, we 
have \r2\ < 1; this means that indeed the remainder term 
that "spoils" the inertial range scaling to first order is 
smaller in the new representation. One can see, by solv- 
ing Eq. (|l^) in terms of the coefficient of the remainder 
term, using -^2^(2^) to express D2n{y) back again and 
comparing with Eq. (12), that one should have the rela- 
tion 



ri 



(24) 



Now comparing Eqs. (15) and ( |23| ) one sees that this 
relation is true. Extended Self Similarity as represented 
by ri very near unity indeed translates into very small 
r2, that is into a small remainder term. In Fig. 2(b) 
we can see r2 as a function of ^2 for n = 2,?i = 4 and 
rt = 6. As n — > 00 the ratio becomes poor, r2 ^ — 1 
but the convergence is slow as we can see from Fig. 2(b) 
and again not uniform in (^2, since one can verify that 
\miQ^^(3r2 = 0. 

The situation for (^2 > 1 is much messier algebraically 
although just the same procedure can be followed in prin- 
ciple. We have already given in Eq. ( p^ ) the near-inertial 
form for D2{y) for (2 > 1; we now substitute this and 
D2n{y) = ^2ny^^" [1 + ^.iv)] [considcrcd as a definition of 
£,{y)] into Eq. (0) to find the equation 



+ (2C2n - 
<2(1 



2C2)(2 



(C2 + 1) 



C2) o 

— r 



(25) 



This equation has an obvious power law solution which 
leads to 



D2n{y) = A2ny^' 



1 



»C2(1 + 2C2) c._2 
(C2 + 1)(1 + 2C2„)^ . 



(for (2 > 1), which gives 

71(2(1 + 2(2) 



ri 



C2n(l + 2C2n) 



In Fig. 3 we can see ri as a function of C2 for n = 2, n 
and n — 6. For very large n we can see that 



lim ri — 



1 



\-_2C2 4 5 
6 ^^2'6^ 



(26) 

(27) 
4 

(28) 



but the convergence is slow as we can see from Fig. 3 
and again not uniform in (2, since one can verify that 
lim^2^o ^1 = 1- Again one can verify that the efficacy of 
ESS is reflected in a smaller remainder term in the ODE 
for S2n using S2 as a variable, as expressed by a ratio r2 
for which here too r2 = 1 — rj~^. 

The success of the method in this case cannot be glibly 
extended to any turbulent system. We see from the anal- 
ysis that the improvement in the scaling behaviour of Sn 
when presented as a function of S2 is not a fundamen- 
tal result but a statement on the similar behaviour of the 
structure functions going into the dissipation range. This 
is measured by the value of the parameter ri of Eq. (^|) . 
We can imagine in principle - although no examples are 
known - a system for which the dissipative scale for Sn 
is smaller than that for 82- In such a system the method 
will not improve the scaling plots. In this passive scalar 
model the dissipative scales of all the structure functions 
are of the same order, as has been shown in [Q. This 
is one of the deep reasons for the success of the method 
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in the present case. It is tempting to conjecture that 
the success of the method in Navier- Stokes turbulence 
is an indication that the dissipative scales of Sn are ei- 
ther of the same order or become larger with n. The full 
understanding of the method in that case must await a 
better analytic understanding of the structure functions 
and their functional form in the near dissipative regime. 
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Figure Legends: 

Fig.l: The dimensionless structre functions D2n{y) 
plotted in log-log plots as functions of y (solid line) and 
S2{y) (dashed-dotted line) respectively. In dashed line 
the exact scaling law is shown. In panels (a) and (b) 
(fi = 1/2 and n = 2 and n = 6 respectively. Panels (c) 
and (d) show results for — 3/2 and the same value so 
of n. One sees an improvement of the scaling law in at 
least 54 orders of magnitude ???. 

Fig. 2: panel (a): The ratio ri of Eq.(15) as a function 
of C2 for n = 2 (solid line) , n = 4 (dashed line) and n = 6 
(dotted-dashed line), panel (b): The ratio r2 of Eq.(25) 
for the same values of n. 

Fig. 3 The ratio ri of Eq.(15) for (2 > 1, and the same 
values of n as in Fig. 2. 
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